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Abstract 

Since several years there exists a question whether the dimensional 
reduction and the usual dimensional regularization give different re- 
sults for the QCD-improved b — > sj and b — » s gluon decay rates. 
Here it is demonstrated explicitly that this is not the case: As long as 
physically meaningful quantities are considered, the results obtained 
with help of both techniques agree. 
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1 Introduction. 



Since several years there exists a question whether the dimensional reduction 
(DRED) and the usual dimensional regularization with fully anticommuting 
7 5 (NDR) give different results for the QCD-improved b —>■ sj and b — > s 
gluon decay rates |J. The purpose of this letter is to show explicitly that this 
is not the case: As long as physically meaningful quantities are considered, 
the results obtained with help of both methods agree. 

The discussion presented here will be based only on the decay b — > sj 
which receives continuous interest at present, because of the recent measure- 
ment of the B — > K*j decay. 

The b — > 57 proceess at the leading order in the Standard Model interac- 
tions is given by the sum of the one- loop diagrams of Fig.l. The QCD con- 
tributions can be diagramatically represented by connecting the quark lines 
of these diagrams with an arbitrary number of gluon lines. The diagrams con- 
structed this way form a power series in the parameter a QCD (M w )ln(M^r / ml) 
~ 0.7 which seems too large to be an expansion parameter. Therefore one 
has to resum all the large logarithms with help of the Operator Product 
Expansion and the Renormalization Group Equations (RGE). 

In order to do this, one introduces the local effective hamiltonianF] 



AG Fl 
~V2 



He ff = 7^ V ts V tb 



Cj(/i)Oj + counterterms 



Li=l 



(1) 



A complete list of the operators Oj is given in each of the papers []3], |5], || [7], 0] . 
As pointed out in ref. f|, the papers fl], |J do not include all the operators 
relevant in the leading-logarithmic approximation. This problem is, however, 
to a large extent unrelated to the problem whether both dimensional methods 
give the same results. 

For the purpose of this paper, it is enough to give explicitly now only 
three of the operators Oj. 

2 = (s L 7 M c L )(c L 7^6 L ) 

7 = -^— 2 m b s L a, u F^b R (2) 

8 = -r^rnbSLd^G^bR 

where F^ u and G^ v are the photonic and gluonic field strenght tensors, re- 
spectively. One finds the coefficients of these operators by requireing equal- 
ity of amplitudes generated by the effective Hamiltonian ([J) and the full 

^or simplicity we neglect the small V* s V u b in our discussion. However, the basic results 
summarized in eqs. (Q), (Jfy and in the Appendix are exactly the same even if V* s V u b is 
not neglected. 
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Standard Model amplitudes (up to 0(1/M^)) at the renormalization scale 
fi = M w . The well-known results are: 



c 7 (M w ) 



c 2 {M w 
3x 3 - 2x 2 



c 8 (M, 



4(x- l) 4 
-3x 2 



w) 



A(x 



ln{x) + 
■ln(x) 4- 



-8a;- 



5a; + 7a; 



24(x - l) 3 
-x 3 + 5a; 2 + 2a; 



,r 



l) 3 



(3) 
(4) 

(5) 



Here x = mf/M^. The above coefficients are regularization- and renorma- 
lization-scheme independent. In order to avoid appearance of large loga- 
rithms in the b — > s'j matrix element of H e ff flU), one evolves the coefficients 
Cj(/x) down to the scale \i = m b , according to the RGE: 



d n 



(6) 



Finally, one finds that the b — *■ 57 matrix element of H e fj (0) evaluated at 
fi = m b (in the leading-logarithmic approximation and for the photon on 
shell) is equal to the tree-level matrix element of fj 



4GV 
72" 



V t * s V tb c e 7 ff (m b )0 7 



(7) 



with 



8 

c/ f {m b ) = V %c 7 {M w ) + - (77M - jyi) c 8 (M w ) + c 2 (M w ) £ a^ 6 * (8) 



i=i 



where 77 = , and a,, 6, are some exact numbers that should be regularization- 

and renormalization-scheme independent. 

The form of the last expression agrees with the final results of any of the 
papers fl], || §, [5], || |7|, ||. However, the particular values of the numbers dj 
disagree for any two of them. The numbers and 6j are given and discussed 
in the Appendix. 

In order to make our discussion as simple as possible, we formally expand 
the last term in in powers of a QCD (Mw)'- 



x ^c»(Mw) ln Mw +Q 



i=i 



2tt 



a QCD (M w )ln 



Mi 



m b 



where 



x = — 2^ aih 

6 i=l 



(9) 



(10) 



2 The s-quark mass is neglected throughout. 
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From the NDR results of refs. @, |, g, g § §, 0, § it follows that X = 
§, , f , =f , ™, *», respectively. The first change (from § to f ) 

is due to the inclusion of the contributions proportional to the down-quark 
charge, that had not been included in the first paper. The second change 
(from ^ to ^p) is due to taking into account certain one-loop matrix el- 
ements which will be described in the next section. The value of has 
been very recently confirmed || with help of the HV scheme where no 
one-loop matrix elements enter. The only known DRED calculation (see ref. 
[|I|]) gives X = i|p In the following two sections, the calculation of X in NDR 
and DRED will be presented and shown to give the same result of We 
restrict ourselves only to the quantity X (which has been the subject of the 
main discussion in the past) in order to avoid considering all the subtleties 
involved in the calculation of the remaining terms in the r.h.s. of eq. (|]) 
where most of the disagreements between the existing NDR calculations are 
located (see the Appendix). 



2 The NDR calculation of X. 

In order to calculate X, we have to trace out all possible leading-logarithmic 
contributions to < sj\H e ff\b > M=mi , that are proportional to C2(Mw)- Simi- 
larly to refs. [|, and differently than in refs. || ||, || |7|, , we will perform 
the calculation without applying the equations of motion to the operators we 
encounter. 

First of all, we have to consider all possible divergent one-loop 1PI dia- 
grams with the O2- vertex, except for those that contain more than one power 
of the QED-coupling e. There are four types of such diagrams, presented on 
Figs. 2a-2d, respectively. 

The diagram in Fig. 2a vanishes for the on-shell photon. This is why there 
is no one-loop mixing between O2 and O7. In effect, the two-loop mixing 
between these two operators becomes important in the leading logarithmic 
approximation. 

All the diagrams in Fig.2b require only counterterms proportional to the 
four-quark operators containing the cl and cl fields. None of such operators 
can have nonvanishing one-loop b — > S7 matrix elements for the on-shell 
photon. Consequently, these diagrams are irrelevant in the calculation of X. 

The diagrams in Fig. 2c generate no divergencies in NDR. 

The divergencies generated by the diagrams in Fig. 2d can be exactly 

3 In this scheme 75 anticommutes with the 4-dimensional 7 M 's, but commutes with the 
remaining ones. It is the only known scheme where problems with traces containing odd 
numbers of 75's do not appear. The treatment of 75 is the only difference between the HV 
and NDR schemes. 
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cancelled by the counterterm (see eq.(|T|)): 



C 2 (fi)Z 2G G (11) 

where 



and 



o =^s L i li (D v G'")b L (12) 
4 

Z 2G = — 7^— — + (terms finite in the limit d — > 4) (13) 

From Z 2G we recover the corresponding element of the anomalous dimension 
matrix 

7 2G = -- + 0(a QCD ) (14) 
and immediately find the relevant term in the solution of the RGE for 

CgO) 

Cain) = c G (M w ) + L 2 (M w )ln^- + 0(a QCD ln^) (15) 

3 fX fX 

We could write all the remaining terms as 0(a QCD ln^ L ) because there is no 
other than 2 operator that has a coefficient of order 1 at fi = My/ and 
mixes with G at zeroth order^ in a QCD . 

Now, let us consider the one-loop on-shell matrix element of G . It is 
given by the diagrams in Fig.3. When a matrix element is considered, we 
have to take into account also the 1PR diagrams. The sum of the diagrams 
in Fig.3 appears to be proportional to the tree-level matrix element of 0-j\ 



Sum of the diagrams \ 4Gf a QCD 2 

in Fig.3 at „ = » t ) = V^^)—- < s 7 \Or\b > 

/ on shell v 

(16) 

Comparing eqs. (0), ©, ©, dUD and ( |16|) one finds the contribution to 



X from the one-loop matrix element of O 



a- 



AX = — — f17l 

^~ matrix element 27 

A contribution to from one-loop matrix elements has been already 
found in ref. |J. As mentioned there, its particular value (that corresponds 
to — ^ above) is correct only in dimensional regularization with fully an- 
ticommuting 75. In any 4-dimensional scheme (or in the HV scheme) the 
on-shell b — > 57 one-loop matrix element of G vanishes. But the contribu- 
tions to C7 from the 2-loop mixings can (or even must) be also different. 
This point has been recently emphasized in ref. M. 



in the applied normalization for O c 
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What remains to be considered is the two-loop mixing O2 — > O7. It is 
described by the diagrams in Fig.4. The calculation of these diagrams is 
described in great detail in ref. ||. The well-known result is 

a QCD 232 

727 = -^TJi (18) 



After solving the RGE of eq. <^)f\ and formally expanding the solution in a QCD 
we immediately get 

232 

^^-^ mixing § X ^ ^ ^ 

Therefore, the final value of X is 

X = AX 4- AX __A,232 = 208 

^/t Z_ivl matrix element ~ T~ mixing c\*-7 oi Ol \ / 

ol ol 

3 The DRED calculation of X. 



The DRED scheme has been introduced by Siegel |Tl[ in order to be able to 
dimensionally regularize supersymmetric theories without actually breaking 
super symmetry. The only difference between NDR and DRED is that in 
the latter scheme all the tensor fields are left 4-dimensional, while the mo- 
menta and coordinates are d-dimensional. The dimension d is assumed to be 
"smaller" than 4, which means that in practical calculations one applies the 
equality^ 

9 ( %9 UP = < (21) 

where g (4 ^ u and g^ u denote the 4-dimensional and the d-dimensional metric 
tensors, respectively. As it is implicit in the above equation, we allow the 
indices of the d-dimensional tensors to acquire values also larger than d, but 
then the corresponding components of these tensors are assumed to vanish. 

For clarity, we will keep the superscript "(4)" for all the 4-dimensional 
tensors that appear in this section. It is also necessary to introduce some 
notation for the diference between the 4-dimensional and the d-dimensional 
tensors. We define: 

G% = G% - G, (22) 

for the gluonic field, and 

1% = - 7, (23) 

for the Dirac matrices. The matrix 75 is taken to be anticommuting with all 
the 7^'s and 7^'s 0. 

We should supply the superscripts "(4)" for all the fields and Dirac ma- 
trices in the definitions of the three operators we started with in eq. @, and 

5 for arbitrary values of all the other mixings 

6 opposite to the HV scheme where g (4 ^g' yp = g 14 ^ is used 
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to the operator in eq. (|I2"D. The "new" operators will be denoted by O 
Of \ Og 4) and respectively. 



(4) 
2 > 



We proceed along the same lines as in the previous section. We have to 
consider the divergent parts of the diagrams in Figs. 2a-2d, but now with 
the O2 -vertex. The diagrams of Figs. 2a and 2b are elliminated from our 
discussion with help of the same arguements as in the NDR case. 

The first real difference between NDR and DRED is in the case of the 
diagrams in Fig. 2c. In the DRED case the sum of the diagrams with one 
gluon and one photon does not vanish, but gives a divergence proportional 
to the operator: 

O x = ^sfa^F^^Hf (24) 
The required renormalization constant is 

4 



'2x 



3(4 - d) 



+ (terms finite in the limit d — *■ 4) (25) 



So, similarly to the case of G discussed in the previous section, we get (cf. 
eqs. dD-dD): 

c x (n) = c x {M w ) - ^c 2 (M w )ln^- + ti{a QCD ln—^-) (26) 

The diagrams in Fig. 2c that contain two gluons are irrelevant in the 
calculation of X, because the one-loop b — > 57 matrix element of an operator 
containing two gluons is of order a^ CD . 

The diagrams in Fig. 2d generate two important counterterms. One of 
them is, of course, the counterterm proportional to with the same renor- 
malization constant as in eq. ([IB]). The other is proportional to the operator 



16vr 

with the renormalization constant 

2 

3(4~ d) 



°y = ih'L^G^ (27) 



2 v = In + (t erms finite in the limit d — > 4) (28) 



So, the coefficient of O y behaves like 

Cyip) = c y (M w ) - \ C2 {M w )ln^- + 0(a QCD ln^) (29) 

The diagrams containing more than one gluon in Fig. 2d may also give 
rise to some other counterterms containing the Gff field. They are, however, 
irrelevant for the same reason as the two-gluon diagrams in Fig. 2c. 
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The appearance of the counterterms involving the that break gauge 
invariance in the (4— d)-dimensional subspace is nothing surprising in DRED. 
They have been observed already by the inventor of DRED [|ll|] . 



Now, let us consider the one-loop on-shell matrix elements of the operators 
Ox-, and O y . The first of them is given by the sum of diagrams in 
Fig. 5. The sum of these diagrams appears to be proportional to the tree 
level matrix element of 7 A) (we ignore the possible terms proportional to the 
(4 — c?)-dimensional photonic field): 



Sum of the diagrams \ 4GV ,a QCD ( 2\ 

V ts V tb c x {m b )^— [--) < s-f\0 7 \b> 



\ in Fig.5 at fi = m b I y/2 2tt V 3 

V / on shell v 

(30) 

Comparing eqs. (0), (§, (g), fl2|) and (^) one finds: 

Q 

A V" _ __ (o-l \ 

matrix element of Ox f\ V / 

y 

The one-loop matrix elements of and O y are described by the same 
diagrams as in the case of G (Fig. 3) - we only have to change the operator 
vertex. In the case of the sum of these diagrams vanishes on-shell, while 
in the case of O y we get exactly the same as in eq. (|16|). The correlation is 
not surprising, because 

G = Of + O y (32) 
Using these results in the same way as before, we get: 



^■■^ matrix element of 0^> ^ matrix element of Oy 

So the sum of the contributions to X from the one-loop matrix elements is 

8 4 20 

^^■X- matrix elements ~ t - ~\~ o*7 (^^0 

y z I z I 



Now, we have to consider the two-loop mixing — > ( 7 4) . The two- 
loop diagrams look exactly the same as in the NDR case (Fig. 4). In the 
one-loop counterterm diagrams of Fig.4 we have to insert both the O^- and 
the O y -counterterms. Finally, we have to take into account also the one-loop 
counterterm diagrams with the counterterms proportional to O x . The latter 
diagrams can be obtained from the ones in Fig.5 just by replacing the square 
(representing the O^-vertex) by a cross (representing the O a .-counterterm). 

The details of the two-loop calculation will not be presented here. It has 



been done with help of the method used in ref. fl3fl , i.e. only the difference 
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between DRED and NDR has been calculated. Then one needs to consider 
only the double-pole parts of the two- loop integrals (given in ref. ||), and 
also the Dirac algebra is relatively simpler. 

The final result for the sought element of the anomalous-dimension matrix 

is: 

a Q cp (124 16\ = a QCD 268 
2vr V 81 9 / 2vr 81 



727 



(35) 



The number comes from the O^-counterterm diagrams. This will be ex- 
actly the contribution from these diagrams to X. It is not an accident that 
it equals to —2 (contribution to X from the one- loop matrix element of O x ). 
This is a common feature for the so-called "evanescent operators" i.e. oper- 
ators vanishing in the limit d — > 4 (see e.g. refs. [14], [15|, |J). 

The number in eq.fl3"5|) comes from the diagrams of Fig.4. It is in 



agreement with the findings of ref. [0, where the O x counterterms have not 
been included. The presence of the O x counterterms is also the reason why 
the tests made in ref. 



121 did not work in the DRED case. 



As in the NDR case, we recover the contribution to X from 7 27 . 



268 
~8T 



(36) 



and we add it to the contribution from the matrix elements, to obtain the 
final result 



v a v a v 20 268 208 

V — A V I A V I _ 

I—\yY matrix elements "I - mixing 0^7 ~ OI 01 

At ol ol 

which is in agreement with the NDR result. 



(37) 



4 Final remarks 

The following table summarizes the results for the quantity X obtained with 
help of the NDR, HV and DRED schemes. 



scheme 


A X 

matrix elements 


A X 

^— mixing 


X 


NDR 


8 

27 


232 
81 


208 

81 


HV 





2U8 

81 


208 

81 


DRED 


20 

27 


268 
81 


208 
81 



The result for AX mixing in the HV scheme has been taken from eq.(25) of 



ref.[ 



The equality of all the three results for the physically meaningful quantity 
X is what one would naturally expect, assuming that all the three schemes 
are the consistent ones. This is also what one could expect by remembering 
the two-loop calculation of ref. [[HJ where all the three schemes were found 
to give the same results for the physically meaningful quantities in the four- 
quark operator case. 

A scheme for extending the Dirac algebra to d-dimensions is consistent 
if it gives a proper limit at d — > 4 and is unique[|. The latter requirement 
follows from the fact that a consistent regularization procedure must give the 
same results for a given diagram independently on whether it is considered 
separately or as a subdiagram, and independently on the order in which the 
subdiagrams are calculated. 

By analyzing how all the three schemes are defined, one can realize that 
all of them are consistent as long as traces with odd numbers of 7s's do not 
appear. If they appear, then only the HV scheme remains consistent. In the 
two other schemes we find that the expressions like 

Tr(7 a 7 At 7^7 p 7 <T 757 Q ) (38) 

in NDR, or 

Tr^^W^l^) (39) 

in DRED give different results dependently on whether the cyclicity of 
the trace is used, or the contracted 7's are commuted all the way through 
what stands between them. 

In the calculation of X presented in this paper, as well as in the calculation 
of ref. [TJJ no traces appeared. This is why no discrepancies between the 



three schemes were observed. 

Things become more complicated in the complete leading-logarithmic cal- 
culation of the b — > 37 rate. In that case one has to use more refined argue- 
ments to show that the (sufficiently careful) NDR calculation is a consistent 
one. Two independent ways of such an argumentation have been given in the 
Appendix A of ref. || . The first was based on the idea of introducing certain 
"evanescent operators" in order to avoid specifying the algebraic properties 
of 75 before arriving at expressions that cannot lead to inconsistencies^. The 
other method was based on the observation, that certain symmetries in the 



7 We do not require that our regularization scheme preserves the symmetries of the 
theory. In the absence of anomalies, all the nonsymmetric terms are local and can be 
removed by proper nonsymmetric counterterms. 

8 This method is a consistent one, but in a general multi-loop calculation it is expected 
to be much more complicated than the usual HV scheme. However, in some particular 
calculations (especially such where no more than two 75 's can appear in a single fcrmionic 
line) it can be much simpler. 
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structure of all the relevant four-quark operators allows to perform the cal- 
culation with no need to calculate any traces. This observation has been also 
independently made in ref . [|ll| , in the context of calculating similar diagrams 
for the next-to-leading effects in the AS = 1 transitions. 

Both these methods of dealing with the dangerous traces can be directly 
used also in the DRED case. This is why one can expect that the sufficiently 
careful DRED calculation of all the numbers etj and 6j from eq.(|) will give 
the same results as in the NDR and HV schemes^. 

The author would like to thank Professor A. Buras for stimulating dis- 
cussions. 



Appendix 

As mentioned below eq.(|J), the numbers Oj and foj present in this equation 
are subject to disagreements between any two of the existing calculations of 
the leading-logarithmic QCD effects in the b — > 57 decay |IJ, |3|, [|, [7], H]. 
In most cases it is due to the fact that not all the relevant dimension-six 
operators are included. The complete basis of the operators (reduced by the 
equations of motion) has been written down already in ref.||. But the effects 
of the operators called[3 O3, O4, O5, Oq have been neglected there. 

The three most recent papers |6|, [7], ||, where the effects of all the rele- 
vant operators are explicitly calculated, disagree on the mixings (0 5 , 6 ) — > 
(0 7 , 0$). The authors of refs.H and M have performed a comparison of 



their calculations [ IT] , and they have found two (and only two) sources of 
the disagreement 

(i) not including the effects of the so-called "evanescent operators" in 
ref.0 

(ii) disagreements in these parts of the first four two-loop diagrams in 
Fig. 5 of ref. || , where the fermion mass comes from the fermionic loop. 

After arriving at this conclusion, we have received the paper ||. The 
results of this paper can be reproduced, if the disagreement (i) is resolved 
in favour of ref. , and the disagreement (ii) - in favour of ref. |?J . However, 
none of the authors of the considered papers is ready to say at the moment 
that his previous results should be corrected. This is why I have decided to 
give here the numbers a, corresponding to each of the three papers [|, [7], §] . 
They are as follows 



9 This is also the expectation of the authors of ref. g who declare to be just performing 
such a complete DRED calculation. 
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The numbering of the operators common for refs. |fj, |) is used here. 
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ref.§: 

Oi = (f|ff , -fff, -0.4286, -0.0714, -0.1991, -0.0453, -0.0215, -0.0990) 
ref.0: 

ai = , -ffff , -0.4286, -0.0714, -0.7415, -0.0003, -0.0580, +0.0323) 

ref.g: 

a, t = (§§§§, -fifi' -0-4286, -0.0714, -0.6494, -0.0380, -0.0186, -0.0057) 

The corresponding numbers hi are the following 

h = (if, g, ^, 0.4086, -0.4230, -0.8994, 0.1456) 

The numbers hi are insensitive to the disagreements between the papers 
IPj • This * s because they are proportional to the eigenvalues of the block- 
triangular anomalous- dimension matrices which disagree with each other only 
in the off-diagonal block. The latter six of the numbers hi are proportional to 
the eigenvalues of the anomalous dimension matrix for the 0\ — Oq oparators 



that has been calculated long ago in ref. fllj|. The first two are given by the 
self-mixing of the Oj and 0$ operators that has been originally found in 
ref.p]. 



Some of the numbers aj and hi are not rational, but they are known 
to arbitrary precision because they come from the diagonalization of the 
leading-order anomalous- dimension matrices that are known exactly. 

As it has been mentioned below eq.@, the sum of all the numbers a« 
always vanishes. This is because all the QCD effects summarized in eq.(|^) 
must vanish for 77 = 1. 

The numbers corresponding to the results of refs.|| and f| look very 
much different. However, the difference between the resulting s is be- 
low 1%. This can be easily undestood, because the differences between the 
anomalous dimension matrices found in these papers are only in the mixings 
O5 ~^ O-j and (0 5 , Oq) — > 0$. The operator 5 acquires only a very small 
coefficient (~ 0.008) during the evolution from Mjy to mj. The operator Oq 
has a larger coefficient, but the effects of the mixing Oq — > Og —>■ Oj tend to 
cancel in the latter step (see the term proportional to 0§ m ec L-(i))- 
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Figure captions 

Fig 1. Diagrams contributing to b — > s'j at the leading order of the SM 
interactions 

Fig 2a-2d. One-loop divergent diagrams with the O2 vertex 
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Fig 3. Diagrams contributing to the one-loop on-shell matrix element 
of G , or O y . The square denotes the insertion of any of these three 
operators. 

Fig 4. Diagrams contributing to the two-loop mixing of 2 with 7 
Fig 5. Diagrams contributing to the one-loop matrix element of O x 
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